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O Abstract 
<N 

^ In this paper we have studied the critical phenomena in higher curvature charged black 

holes in the anti-de Sitter (AdS) space-time. As an example we have considered the third 
order Lovelock-Born-Infeld black holes in AdS space-time. We have analytically derived the 
thermodynamic quantities of the system. Our analysis revealed the onset of a higher order 
phase transition in the black hole leading to an infinite discontinuity in the specific heat at 
constant charge at the critical points. Our entire analysis is based on the canonical framework 
where we have fixed the charge of the black hole. In an attempt to study the behavior of 
the thermodynamic quantities near the critical points we have derived the critical exponents 
of the system explicitly. Although the values of the critical points have been determined 

5— ( numerically, the critical exponents are calculated analytically. Our results fit well with the 

thermodynamic scaling laws. The scaling hypothesis is also seen to be consistent with these 
scaling laws. We find that all types of AdS black holes, studied so far, indeed belong to the 
I same universality class. Moreover these results are consistent with the mean field theory 

approximation. We have derived the suggestive values of the other two critical exponents 
' J associated with the correlation function and correlation length on the critical surface. 

1 Introduction 

in 

o 

String theory and brane world cosmology underscored the necessity of considering space-time 
as multidimensional. This motivates us to study gravity theories in higher dimensions (i.e. 
greater than four)[lj. The effect of string theory on gravity may be understood by considering 
a low energy effective action that describes the classical gravity[2j. This effective action must 
include combinations of the higher curvature terms and are found to be ghost free [3]. In an 
attempt to obtain the most general tensor that satisfies the properties of Einstein's tensor in 
higher dimensions, Lovelock proposed an effective action that contains higher curvature terms[l]. 
Moreover, the field equations derived from this action consists of only second derivatives of the 
metric and hence are free of ghosts [S]. The study of black holes in higher dimensions has found 
renewed attention because of the enriched physics associated with them[6j-[8j. 

Thermodynamics of black holes has been a principal topic of research due to the discovery 
of the laws of black hole mechanics [9]. The study of thermodynamic properties of black holes 
in anti-de Sitter space-time has got renewed attention due to the discovery of phase transition 
in Schwarzschild-AdS black holes[10j. Since then a wide variety of research has been done to 
study the phase transition in black holes[71 [8], [H]-[22]. Now a days, the study of thermody- 
namics of black holes in AdS space-time is very much important in the context of AdS/CFT 
duality. The thermodynamics of AdS black holes may en-light the underlying phase structure 
and thermodynamic properties of CFTs[33]. The study of thermodynamic phenomena in black 
holes requires an analogy between the variables in ordinary thermodynamics and those in black 
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hole mechanics. Recently, R. Banerjee et. al. have developed a method [T8] IT9"] . based on 
Ehrenfest's scheme [23] of standard thermodynamics, to study the phase transition phenomena 
in black holes. In this approach one can actually determine the order of phase transition once 
the relevant thermodynamic variables are identified for the black holes. This method has been 
successfully applied in the four dimensional black holes in AdS space-time |19j- [22] as well as 
in higher dimensional AdS black holes [7]. While most of the black holes undergo second order 
phase transitions characterized by an infinite discontinuity in the specific heat at the critical 
points (i.e. points of divergence), one exception is observed in the charged Hofava-Lifshitz black 
holes in AdS space [25]. Here we observe both first and second order phase transitions. 

The behavior of the thermodynamic variables near the critical points can be studied by 
means of a set of indices known as static critical exponents [261 127] . These exponents are to 
a large extent universal, independent of the spatial dimensionality of the system and obey 
thermodynamic scaling laws|26} 127] . The critical phenomena has been studied extensively in 
familiar physical systems like the Ising model (two and three dimensions), magnetic systems, 
elementary particles, hydrodynamic systems etc.. An attempt to study the critical phenomena 
in black holes was commenced in the last twenty years [28 -[43J. Despite all these attempts, a 
systematic study of critical phenomena in black holes was still lacking. This problem has been 
overcome very recently [HI HI]- In these works the critical phenomena have been studied in (3+1) 
as well as higher dimensional AdS black holes. Also the critical exponents of the black holes 
have been determined by explicit analytic calculations. 

All the research mentioned above was confined to Einstein gravity. It would then be inter- 
esting to study gravity theories in which action involves higher curvature terms. Among the 
higher curvature black holes, Gauss-Bonnet and Lovelock black holes will be suitable candidates 
to study. The thermodynamic properties and phase transitions have been studied in the Gauss- 
Bonnet AdS (GB-AdS) black holes [15]- [19] . Also, the critical phenomena in the GB-AdS black 
holes was studied[50j. On the other hand, Lovelock gravity coupled to the Maxwell field was 
investigated in |51| [52] . The thermodynamic properties of the third order Lovelock-Born-Infeld- 
AdS (LBI-AdS) black holes were studied in [51], [531 E3] - But the study of phase transition and 
critical phenomena have not been done yet in these black holes. 

In this paper we have studied the critical phenomena in the third order LBI-AdS black 
holes in seven dimensions. We also have given a qualitative discussion about the possibility of 
higher order phase transition in this type of black holes. We have determined the static critical 
exponents for these black holes and showed that these exponents obey the static scaling laws. 
We have also checked the static scaling hypothesis and calculated the scaling parameters. We 
see that these critical exponents take the mean field values. From our study of the critical 
phenomena we may infer that the third order LBI-AdS black holes yield results consistent with 
the mean field theory approximation. There are some distinct features in the higher curvature 
gravity theory. For example, the usual area law valid in Einstein gravity does not hold in these 
gravity theories. However, the critical exponents are identical with those in the usual Einstein 
gravity. This result shows that the AdS black holes belong to the same universality class. 
We have also determined the critical exponents associated with the correlation function and 
correlation length. However the values of these exponents are more suggestive than definitive. 
As a final remark, we have given a qualitative argument for the determination of these last two 
exponents. 

Let us mention about the outline of the paper. In section 2 we discuss about the thermo- 
dynamical variables of the seven dimensional third order LBI-AdS black holes. In section 3 we 
analyze the phase transition and stability of these black holes. The critical exponents, scaling 
laws and scaling hypothesis are discussed in section 4. Finally, we have drawn our conclusions 
in section 5. 



2 



2 Thermodynamic variables of higher curvature charged AdS 
black holes 



The effective action in the Lovelock gravity in (n + 1) dimensions can be written a^jl"] 



i=0 



where a, is an arbitrary constant and Li is the Euler density of a 2i dimensional manifold. In 
(n + 1) dimensions all terms for which i > [(n + l)/2] are equal to zero, the term i = (n + l)/2 
is a topological term and terms for which i < [(n + l)/2] contribute to the field equations. 
Since we are studying third order Lovelock gravity in the presence of Born-Infeld nonlinear 
electrodynamics |55j the effective action (1) may be written as 

1 = lib / dn+lx ^9( a o£o + a 1 £ 1 + a 2 C 2 + a 3 £ 3 + L(F)) 

= \kr / dn+lx ^9 (-2A + 11 + a 2 C 2 + a 3 j0 3 + L(F)) (2) 

where A is the cosmological constant given by — n(n — l)/2l 2 , I being the string length, a 2 and 
«3 are the second and third order Lovelock coefficients, C\ = 7Z is the usual Einstein-Hilbert 
Lagrangian, £-2 = R^&R^ 5 - \R UV W + TZ 2 is the Gauss-Bonnet La grangian and 



— 2R iiVaK R RP T -\-RR^ u R aK R pT + 2AR tMVaK R R p 4- 
3nR ^™ R + 24R» uaK R au R K u + 16R^ u RuaR a * - VITIR^ R, v + V? 



is the third order Lovelock Lagrangian, L(F) is the Born-Infeld Lagrangian given by 

l w = 462 ( x - f+w) ^ 

where F^ u = d^A u — d v A^^ F 2 = F^F^" and b is the Born-Infeld parameter. In the limit 
b — > oo we recover the standard Maxwell form L(F) = —F 2 . 

The solution of the third order Lovelock-Born-Infeld anti de-Sitter black hole (LBI-AdS) in 
(n + l)-dimensions can be written as[ 



where dO? n _ x 



ds 2 = -f{r)dt 2 + J^dr 2 + r 2 dtit n _ x (4) 

del + E?=2 n}=i sin%dQ\ k = 1 

d9\ + sinh 2 Q x dB\ + sinh 2 e 1 YZ=s YYj^l sin 2 0jd9 2 k = -1 
I YZi dtf ' ' k = 

and A; determines the structure of the black hole horizorj£j 

Since we are considering third order LBI-AdS black hole, we shall restrict ourselves in the 
seven dimensional space-time (n + 1 = 7). Moreover, the Lagrangian of (2) is the most general 
Lagrangian in seven space-time dimensions that produces the second order field equations [5 lj. 



1 Here we have taken the gravitational constant G = 1. 
2 k — +l(spherical), — l(hyperbolic), 0(planar) 
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One can show that in the special case «3 = 2a\ = y^, the metric function f(r) of (4) may be 
expressed as 



where 



X(r) = 1 + 



f(r) = k 



3am 2ab 2 



a 



1 - x(r) 



In (6) 



Wq 2 
52 r io 



and %(??) is the hyper geometric function given by|56j 



(5) 
(6) 

(7) 
(8) 



q and m are constants which are related to the charge Q and the ADM mass M of of the 
black hole through the equations 



Q 



and 



K__i./(n-l)(n-2) 

47T 



M= -5=^(n-l)m 
16vr V y 



10vr^ 



5^ 
16" 



-77? 



(9) 
(10) 



respectively. Here V n -i = r(n/2) ^ s * ne v °l ume °f the ( n — 1) sphere. 

The ADM mass (M) of the third order LBI-AdS black hole may be defined as f(r+) = 0, 
which yields, 



M 



5tt z 



a 2 4 2 26 2 r6 
3 + + 15 



A 



1-^T+^ + 



20Q 2 
2b 2 1 6W 



(11) 



The electrostatic potential difference between the black hole horizon and the infinity may be 
defined as! 



where 



2(n - 2) 



V-i 



16Q 2 
6 2 7T 4 rI° ' 



(12) 
(13) 



Here we have used (9). 

Using the metric function (5) one may find the Hawking temperature for the third order 
LBI-AdS black hole as, 

T = J_fdf(r] 
4ir 



Or 



Wkrl + hkar\ + 2b 2 r% (1 - ^1 + 77+) - Ar% 
10irr + (r+ + ka) 2 



(14) 



In the limit a — > the corresponding expression for the Hawking temperature of the Born- 
Infeld AdS (BI-AdS) black hole can be recovered as, 



Tsi-AdS 



1 

Air 



4k Ab 2 r A 
r + 5 



1 



'1 + 



Q 2 

b 2 r\_ 



2Ar_, 



(15) 
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From the first law of black hole thermodynamics dM = TdS + &dQ, we may obtain the 
entropy of the third order LBI-AdS black hole as, 





vr 3 / K Wkarl 

T (r: 



5 +^7++5fc 2 a 2 r + . (16) 



From (16) we find that the entropy is not proportional to the one-fourth of the horizon area 
as in the case of the black holes in the Einstein gravity. This is the most striking feature of 
the higher curvature gravity, the usual area law of black hole [57J breaks down in these gravity 
theories. If we take the limit a — > 0, we can recover the usual area law of black hole entropy in 
the BI-AdS black hole [8] 

S=^r\. (17) 

In our study of critical phenomena we will be mainly concerned about the spherically sym- 
metric space-time. In this regard we will always take the value of k to be +1. Substituting 
k = 1 in (14) and (16) we finally obtain the expressions for the Hawking temperature and the 
entropy of the third order LBI-AdS black hole as, 



lOrl + 5arl + 2b 2 r% 1 - Jl + ) - Ar, 6 



and 



T = ) 9 v - (18) 

10irr+(r$. + a)* 



3 Phase transition and stability of the third order LBI-AdS 
black hole 

In this section we aim to discuss the nature of phase transition and the stability in the third 
order LBI-AdS black hole. The nature of phase transition of a thermodynamic system is usually 
determined by inspecting into the discontinuity in the relevant thermodynamic quantities at the 
critical point |24j. While the first order phase transition is characterized by the discontinuity in 
the temperature of the system, any discontinuity in the specific heat gives rise to a higher order 
phase transition. Since black holes are considered as thermodynamic objects one may also expect 
phase transition in these systems. Earlier research works revealed that phase transitions indeed 
occur in black holes [7J [8J, |10]-|22j. In order to determine the order of phase transition in black 
holes we generally use the Ehrenfest's equations for black holes where the usual thermodynamic 
variables are replaced by those for the black holes [IB]. We can further explore the stability 
of the third order LBI-AdS black hole by analyzing the specific heat of the black hole. The 
negative value of specific heat gives rise to the unstable configuration of the black hole whereas 
the positive value of the specific heat makes the black hole stable. We can further clarify the 
stability analysis of the black hole by discussing the nature of the free energy of the black hole. 

Let us now plot the Hawking temperature (T) of the third order LBI-AdS black hole against 
the radius of the outer horizon (r + ) of the black hole. 

We observe from the above figures that, there is a 'hump' and a 'dip' in the T — r + graph. 
Also it is very interesting to note that, the graphs are continuous in r + . Since the entropy (S) of 
the black hole increases with the radius of the outer horizon (r+ ) , we may argue that the entropy 
of the third order LBI-AdS black hole is also a continuous function. This rules out the possibility 
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Figure 1: Plot of Hawking temperature (T) Figure 2: Plot of Hawking temperature (T) 
against horizon radius (r+), for a = 0.5, Q = against horizon radius (r + ), for a = 1.0, Q = 
0.50 and b = 10. 0.50 and b = 10. 



of first order phase transition. In our subsequent analysis we will see that the specific heat at 
constant charge (Cq) shows an infinite discontinuity at the critical points which correspond to 
the 'hump' and 'dip' of the T — r + plots respectively. Without going into the algebraic details, we 
can argue, from our previous work [22 J, that this corresponds to a continuous higher order phase 
transition. In order to study the critical phenomena in the canonical ensemble framework, we 
need to calculate the specific heat at constant charge (Cq). Using the standard thermodynamic 
definition we may write the specific heat at constant charge of the third order LBI-AdS black 
hole as, 



Cr 



T 



df 
,(dS/dr+) 



Q 



where 



Af(r + ,Q) = Wrl(r 2 + + a) 3 



(dT/dr + ) c 
V(r+,Q) 



10r 2 + 5a + 2b 2 ri 



'1 + 



16Q 2 
6 2 vr 4 ri° 



(20) 



(21) 



and 



V(r+,Q) = 128Q 2 + [l5vr 4 r^a + h^r\ a 2 - Avr 4 r^° - 5vr 4 r^(2 + aA)} 
l 1 + J 6 ^! _ (2& 2 vr 4 rf + lO&Vr* a) (l - Jl+ 16Q ' ^ 



6 2 vr 4 rl° 



(22) 



where we have used (18) and (19). 

The analytic solution of the denominator T>(r+,Q) of (20) is very much problematic for 
arbitrary choices of the parameters a, b and Q. In order to see whether there is any bound in 
the values of b and Q we will consider the extremal third order LBI-AdS black hole. In this case 
both f(r) and J- vanish at the degenerate horizon r e . The above two conditions for extremality 
results the following equation, 



Wrt + 5ar 2 + 2b 2 rf 1 



1 + 



16Q 2 
6 2 7r 4 ri° 



Art 



(23) 
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Table la: Roots of equation (23) for a = 0.5 and I = 10 



Q 


b 


r 2 
r el 


r e2 


r' 1 

'e3,e4 


r 1 
r e5 


15 


0.6 


-66.4157 


- 


- 


+0.632734 


8 


0.2 


-66.4157 


- 


- 


+0.121590 


5 


0.5 


-66.4157 






+0.200244 


0.8 


20 


-66.4157 


-0.408088 


-0.05881+^0.304681 


+0.268514 


0.3 


15 


-66.4157 


-0.304225 


-0.0517567±i0.175488 


+0.145685 


0.5 


10 


-66.4157 


-0.349764 


-0.0593449±i0. 240144 


+0.192519 


0.5 


1 


-66.4157 






+0.0221585 


0.5 


0.5 


-66.4157 






+0.00625335 


0.05 


0.05 


-66.4157 






+6.57019xl0- ? 



Table lb: Roots of equation (23) for a = 1.0 and I = 10 



Q 


b 




r e2 


r e3,e4 


' e5 


15 


0.6 


-66.1628 






+0.50473 


8 


0.2 


-66.1628 






+0.0546589 


5 


0.5 


-66.1628 






+0.110054 


0.8 


20 


-66.1629 


-0.563564 


-0.0913395±i0.266887 


+0.236186 


0.3 


15 


-66.1629 


-0.514815 


-0.0584676+^0.14609 


+0.116834 


0.5 


10 


-66.1629 


-0.531635 


-0.0760937±i0.210702 


+0.155024 


0.5 


1 


-66.1629 


-0.542417 




+0.00640486 


0.5 


0.5 


-66.1629 






+0.00163189 


0.05 


0.05 


-66.1629 






+1.64256xl0~ 7 



In Table la and lb we give the numerical solution of (23) for different choices of the values of 
the parameters b and Q for fixed values of a. From these analysis we observe that for arbitrary 
choices of the parameters b and Q we always get atleast one real positive root of (23). This 
implies that there exists a smooth extremal limit for arbitrary b and Q and there is no bound 
on the parameter space for a particular value of a. This result is indeed expected since it was 
observed that, although, there is a bound on the parameters b and Q for four dimensional (n = 3) 
Born-Infeld AdS black holes, the bound is removed if we consider these black holes in higher 
dimensions (n > 3) [8| 135]. 

We shall now find out the roots of (22) for different values of the parameters b, Q and a. 
The numerical results are given in the Tables 2a and 2b below. For convenience we shall write 
the real roots of (22) only. From our analysis it is worth mentioning that the specific heat 
at constant charge of the third order LBI-AdS black hole always possesses simple poles. Also, 
there are two real positive roots (ri, r^) of the denominator of Cq for different values of the 
parameters b, Q and a. Moreover, these two roots (r\ and r%) correspond exactly to the 'hump' 
and the 'dip' in the T — r + plot. Another interesting point can be observed from the Cq — r + 
plots that Cq shows an infinite discontinuity at the points n and r^, which may be identified as 
the critical points of phase transition in the third order LBI-AdS black hole. From the nature 
of the T — r+ and Cq — r+ plots it is indeed obvious that the phase transition in the third order 
LBI-AdS black hole is not a first order phase transition, rather, the discontinuity of Cq at the 
critical points {r\ and r2) reveals that it is a continuous higher order (i.e. greater than first) 
phase transition. 

We can further analyze the stability of the third order LBI-AdS black hole once we study 
the behavior of Cq at the critical points. The Cq — r + plot shows that there are indeed three 
phases of the black hole. These phases can be classified as - Phase I (0 < r + < r%), Phase 
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Table 2a: Real roots of equation (22) for a = 0.5 and I = 10 



Q 


b 


n 


7*2 


^3 


r 4 


15 


0.6 


1.59399 


7.96087 


-7.96087 


-1.59399 


8 


0.2 


1.18048 


7.96088 


-7.96088 


-1.18048 


5 


0.5 


1.25190 


7.96088 


-7.96088 


-1.25190 


0.8 


20 


1.01695 


7.96088 


-7.96088 


-1.01695 


0.3 


15 


0.975037 


7.96088 


-7.96088 


-0.975037 


0.5 


10 


0.989576 


7.96088 


-7.96088 


-0.989576 


0.5 


1 


0.989049 


7.96088 


-7.96088 


-0.989049 


0.5 


0.5 


0.987609 


7.96088 


-7.96088 


-0.987609 


0.05 


0.05 


0.965701 


7.96088 


-7.96088 


-0.965701 



Table 2b: Real roots of equation (22) for a = 1.0 and / = 10 



Q 


b 


n 


T2 


^3 


r 4 


15 


0.6 


1.76824 


7.74534 


-7.74534 


-1.76824 


8 


0.2 


1.53178 


7.74535 


-7.74535 


-1.53178 


5 


0.5 


1.51668 


7.74535 


-7.74535 


-1.51668 


0.8 


20 


1.40553 


7.74535 


-7.74535 


-1.40553 


0.3 


15 


1.40071 


7.74535 


-7.74535 


-1.40071 


0.5 


10 


1.40214 


7.74535 


-7.74535 


-1.40214 


0.5 


1 


1.40214 


7.74535 


-7.74535 


-1.40214 


0.5 


0.5 


1.40213 


7.74535 


-7.74535 


-1.40213 


0.05 


0.05 


1.39992 


7.74535 


-7.74535 


-1.39992 



II (n < r + < r<i) and Phase III (r+ > r-i). Since the higher mass black hole possesses larger 
entropy /horizon radius, there is a phase transition at t\ from smaller mass black hole (Phase 
I) to intermediate (higher mass) black hole (Phase II). The critical point T2 corresponds to a 
phase transition from an intermediate (higher mass) black hole (Phase II) to a larger mass black 
hole (Phase III). Moreover, from the Cq — r + plots we note that the specific heat Cq is positive 
for Phase I and Phase III whereas it is negative for Phase II. Therefore Phase I and Phase 

III correspond to thermodynamically stable phase (Cq > 0) whereas Phase II corresponds to 
thermodynamically unstable phase (Cq < 0). 

We can further extend our stability analysis by considering the free energy of the third order 
LBI-AdS black hole. The free energy plays an important role in the theory of phase transition 
and critical phenomena. We may define the free energy of the third order LBI-AdS black hole 
as, 

F(r + ,Q) = M(r + ,Q)-TS (24) 

where M(r+, Q), T and S are the ADM mass, the Hawking temperature and the entropy of the 
black hole respectively. 

Using (11), (18) and (19) we can write (24) as 



5vH 
~16 



a 



+ r , + ar, + 



2b 2 r 



2^6 



15 



1 + 



16Q 2 



A 



6% 4 r|° 2b 2 6 2 tt 4 



7r 2 (rl + ^± + 5a 2 r + ) 



10r^ + har\ + 2b 2 rl I 1 




40r+(r + + a) 2 

We plot the free energy (T) of the black hole with the radius of the outer horizon r + in the 
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Figure 3: Plot of specific heat (Cq) against hori- Figure 4: Plot of specific heat (Cq) against hori- 
zon radius (r + ), for a = 0.5, Q = 15 and zon radius (r + ), for a = 0.5, Q = 15 and 
b = 0.60 at the critical point r\. b = 0.60 at the critical point r-i- 




Figure 5: Plot of specific heat (Cq) against hori- Figure 6: Plot of specific heat (Cq) against hori- 
zon radius (r + ), for a = 0.5, Q = 0.50 and zon radius (r + ), for a = 0.5, Q = 0.50 and 
b = 10 at the critical point r\. b = 10 at the critical point r2- 



following figures (Fig. 11-14). 

Let us now analyze the J- —r + plots. The free energy (F) has a minima T = T m at r + = r m . 
This point of minimum free energy is exactly the same as the first critical point r + = r\ where 
the black hole shifts from a stable to an unstable phase. On the other hand T has a maxima 
T = J-q at r+ = ro . The point at which T reaches its maximum value, is identical with the 
second critical point r + = r<i where the black hole changes from unstable to stable phase. We 
can further divide the T — r+ plot into three distinct regions. In the first region r[ < r+ < r m 
the negative free energy decreases until it reaches the minimum value (J- m ) at r + = r m . This 
region corresponds to the stable phase (Phase I: Cq > 0) of the black hole. The free energy 
changes its slope at r + = r m and continues to increase in the second region r m < r + < ro 
approaching towards the maximum value (J~o) at r + = ro. This region corresponds to the Phase 
II of the Cq — r + plot where the black hole becomes unstable {Cq < 0). The free energy changes 
its slope once again at r + = ro and decreases to zero at r + = r' 2 and finally becomes negative 
for r + > r' 2 . This region of the J- — r + plot corresponds to the Phase III of the Cq — r_|_ plot 
where the black hole finally becomes stable (Cq > 0). 
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Figure 7: Plot of specific heat (Cq) against hori- 
zon radius (r + ), for a = 1.0, Q = 5 and b = 0.5 
at the critical point r\ . 
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Figure 9: Plot of specific heat (Cq) against hori- 
zon radius (r + ), for a = 1.0, Q = 0.05 and 
b = 0.05 at the critical point r\. 
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Figure 8: Plot of specific heat (Cq) against hori- 
zon radius (r + ), for a = 1.0, Q = 5 and b = 0.5 
at the critical point r2 ■ 




Figure 10: Plot of specific heat (Cq) against 
horizon radius (r + ), for a = 1.0, Q = 0.05 and 
b = 0.05 at the critical point ri- 



4 Critical exponents and scaling hypothesis 

In thermodynamics, the theory of phase transition plays a crucial role to understand the behavior 
of a thermodynamic system. The behavior of thermodynamic quantities near the critical point (s) 
of phase transition gives a considerable amount of information about the system. The behavior 
of a thermodynamic system near the critical point(s) is usually studied by means of a set of 
indices known as the critical exponents [261 EZ]- These are generally denoted by a set of Greek 
letters: a, f3, 7, 5, <p, if), n, v. The critical exponents describe the nature of singularities in 
various measurable thermodynamic quantities near the critical point(s). 

In this section we aim to determine the first six static critical exponents (a, (3, 7, 5, (f), ip). 
We shall then discuss about the static scaling laws and static scaling hypothesis. We shall 
determine the other two critical exponents (y and rj) from two additional scaling laws. 

Critical exponent a: 

In order to determine the critical exponent a which is associated with the singularity of Cq near 
the critical points rj (i = 1, 2), we choose a point in the infinitesimal neighborhood of as, 

r + = r<(l + A), i = l, 2 (26) 

where |A| << 1. Let us denote the temperature at the critical point by T(rj) and define the 
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Figure 11: Plot of free energy (J 7 ) against hori- Figure 12: Plot of free energy(T) against hori- 
zon radius (r + ), for a = 0.50, Q = 0.50 and zon radius (r + ), for a = 0.50, Q = 0.50 and 
6 = 10. 6=10. 





Figure 13: Plot of free energy(T) against hori- Figure 14: Plot of free energy(T) against hori- 
zon radius (r + ), for a = 1.0, Q = 5 and b = 0.5. zon radius (r + ), for a = 1.0, Q = 5 and b = 0.5. 



quantity 



T(r + )-T{ n ) 

Tin) 



(27) 



such that |e| << 1. 

We now Taylor expand T(r+) in the neighborhood of rj keeping the charge constant (Q 
Q c ), which yields, 



T{r + ) = T{ n ) + 



dT 

dr+ 



(r+-ri) + 



1 



J r+=n 



d 2 T 
dr\ 



( r+ -n) 2 + (28) 



r+=ri 



Since the divergence of Cq results from the vanishing of ) at the critical point r, ((20)), 
we may write (28) as, 



T(r + ) = T(n) + - 



d 2 r 



Q=Qc 



(r + - rif 



(29) 



r+=r; 



where we have neglected the higher order terms in (28). 



11 



Using (26) and (27) we can finally write (29) as, 

e V2 



A = 



,1/2 



where 



2T{ n ) 



8 2 T 
drl 



(30) 



(31) 



r+=n 



The detailed expression of Tj is very much cumbersome and we shall not write it down for 
the present work. 

If we examine the T — r+ plots (Fig. 1 and 2) we observe that, near the critical point r + = r\ 
(which corresponds to the 'hump') T(r+) < T(n) so that e < and on the contrary, near the 
critical point r + = r 2 (which corresponds to the 'dip') T(r + ) > T(r2) implying e > 0. 

Substituting (26) into (20) we can write the singular part of Cq as, 



C Q = 



N'(r u Q c ) 
A-V>(r h Q c ) 



(32) 



where M'(ri, Q c ) is the value of the numerator of Cq ((21)) at the critical point r + = rj and 
critical charge Q = Q c . The expression for V(ri, Q c ) is given by, 



v'{ n , Q c ) = vi(n, Q c ) + v' 2 (n, Q c ) + v^n, Q c ) 



(33) 



where 



V' 1 (r i ,Q c ) = 107r 4 rtJl + 



16Q 2 



b 2 ir 4 r}° 



[{2a 2 + 26 2 rf + 8b 2 rfa) - (Arf + 9ar 2 + 4(2 + Aa)rf )] 



V' 2 {n,Q c ) 



80Q 2 
b 2 r 2 



■D' 3 (n,Q c ) = 20b 2 iT 4 rf 



15a 5a z 

0~ ~A~ 



1 + 



Ar 2 - 5(2 + Aa) 



b 2 ir 4 r}° 



+ 4a 1 + 



10Q 2 

6 2 7r 4 r 10 



It is to be noted that while expanding the denominator of Cq we have retained the terms 
which are linear in A and all other higher order terms of A have been neglected. 

Using (32) we may summarize the critical behavior of Cq near the critical points [r\ and r 2 ) 
as follows: 





A 




: (-6)V2 ; 


[ 






L (+6)V2 J 



e < 



Ti=r\ 



ri=r 2 



€ > 



where 



(34) 



(35) 



V>(n,Q c ) 

We can combine the r.h.s of (34) into a single expression, which describes the singular nature 
of Cq near the critical point rj, yielding 



Cq 



11/2 



1/2 



|r-Tj|V2- 



(36) 
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where we have used(27). Here we have abbreviated T and Tj for T(r + ) and T(ri) respectively. 
We can now compare (36) with the standard form 



C Q ~\T-T t \~ a 



(37) 



which gives a = -. 
Critical exponent j3 : 

The critical exponent (3 is related to the electric potential at infinity (<&) by the relation 

*(r + )-*(r i )~|T-r i |^ (38) 

where the charge (Q) is kept constant. 

Near the critical point r+ = the Taylor expansion of 3>(r+) yields, 



$(r+) = $(n) + 



dr+ 



Q=Qc 



(r+ - ri) + 



(39) 



r+=ri 



Neglecting the higher order terms and using (12) and (30) we can rewrite (39) as, 



$(r+) - *(ri) 



4Q C 



\ 



7r 2 r 4 r V2 r 1 / 2 A 4- 16Q " 



IT - T,-| 1/2 



(40) 



Comparing (40) with (38) we finally obtain /3 = -. 
Critical exponent 7: 

We shall now determine the critical exponent 7 which is associated with the singularity of 
the inverse of the isothermal compressibility (k,^, 1 ) at constant charge Q = Q c near the critical 
point r + = ^ as, 

K T l ~\T-Ti\^ (41) 
In order to calculate /-c^ 1 we use the standard thermodynamic definition 



Q 



-Q 



dQ J T 

(d$\ fdT\ 



\dT) {dQ) 



(42) 



where in the last line of (42) we have used the identity 

cM>\ /dT\ (dQ\ 



dTj Q \d<z>J T 

Using (12) and (14) we can write (42) as, 

-1 fi(r+,Q) 



(43) 



(44) 
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where V(r + ,Q) is the denominator identically equal to (22) (the denominator of Cq) and 
the expression for £l(r+, Q) may be written as, 



n( r+ , Q) = -9—[i28Q 2 + (157rVja + 5vr 4 r 4 a 2 - Avr 4 ^ - 5vr 4 r^(2 + aA))i/l + (16Q 2 /^M°) 
-(26 2 7r 4 ri° + lO&Vr® a) (l - yjl + (16Q 2 /& 2 7r 4 r!°)^ W(2/5, 1/2, 7/5, -16Q 2 /b\ A r™)} - E(r+, Q)(45) 



where 



E(r+,Q) = -^[2&Vr 4 (r 2 + 5a) + (16Q 2 /6 2 7r 4 ri°) + 



5 

vr 2 (-15r 2 a - 5a 2 + r%(A - 2b 2 ) + 5r 4 (2 + Aa - 26 2 a))] 

From (44) we observe that k^ 1 possesses simple poles. Moreover k^ 1 and Cq exhibit common 
singularities. 

We are now interested in the behavior of k^ 1 near the critical point r + = r«. In order to 
do so we substitute (26) into (44). The resulting equation for the singular part of k t 1 may be 
written as, 

x _ Df(r it Q e ) 
Kt ~ A.V'(r t ,Q c y (46) 

In (46), Q'(ri,Q c ) is the value of the numerator of k^ 1 ((45)) at the critical point r + = 
and critical charge Q = Q c . Whereas V(ri,Q c ) was identified earlier ((33)). 

Substituting (30) in (46) we may express the singular nature of k^ 1 near the critical points 
(ri and r%) as, 



IX rp 



where 





Bi 1 




L(- e )V»J 


< 






L(+e)VaJ 








V(r 



n=ri 



€ < 

e > 



(47) 



(48) 



Combining the r.h.s of (47) into a single expression as before, we can express the singular 
behavior of k^ 1 near the critical point r« as, 



:|l/2 



1/2 



IT-TA 1 / 2 ' 



(49) 



Comparing (49) with (41) we find 7 = -. 
Critical exponent 5: 

Let us now calculate the critical exponent 5 which is associated with the electrostatic poten- 
tial ($) for the fixed value T = T% of temperature. The relation can be written as, 



<D(r+) - $( ri ) ~ \Q - Q 



(50) 
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in this relation Qi is the value of charge (Q) at the critical point rj. In order to obtain S we first 
Taylor expand Q(r+) around the critical point r + = r; L . This yields, 



Q(r+) = Q(r t ) + 



dQ 
<9r_i_ 



T=T t 



(r+ - n) + 



1 



(PQ 

d r + / T=Ti 



( r+ - n ) 2 + (51) 



r+=n 



Neglecting the higher order terms we can write (51) as, 



Q(r+) - Q{n 



drl 



(r+ - nY 



(52) 



r+=ri 



Here we have used the standard thermodynamic identity 



[(f) ' 




( dr+\ 


\dr+J T _ 


r+=ri 





r+=ri 



\dQj 



(53) 



r+=ri 



and considered the fact that at the critical point r + = rj, (^~§^~^j vanishes. 
Let us now define a quantity 



T = 



Q(r+) -Qi Q-Qi 



Qr 



Qi 



(54) 



where |T| << 1. Here we denote Q(r+) and Q(ri) by Q and Qi respectively. 
Using (26) and (54) from (52) we obtain 



A = 



T l/2 



1/2 



2Q l 



1/2 



(55) 



where 



*i = 



(56) 



r + =r i 



The expression for is very much cumbersome and we shall not write it in this paper. 
We shall now consider the functional relation 



$ = $(r+,Q) 



from which we may write, 



TJ 



r+=n 



dr+ 



+ 



r+=n 



\dr+) T _ 



r+=n 



dQ 



(57) 



(58) 



r+=r t 



Using (53) we can rewrite (58) as, 

<9$ 



T=Ti 



r+=ri 



dr+ 



Q=Qc 



(59) 



r+=n 



Now the Taylor expansion of $ at constant temperature around r+ = n yields, 

<93> 



$(r+) = $(fi) + 



dr i 



(r+ - n 



(60) 



r+=r-j 
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where we have neglected all the higher order terms. 

Finally using (55), (59) and (12) we may write (60) as, 



$(r+) - $(ri) 



-4Q C 



\ 



7rV 5 Jl + 



16Q2 



2 \ 2 , ,1 



(61) 



& 2 7r 4 r4° / 



Comparing (50) and (61) we find that 5 = 2. 
Critical exponent (f>: 

The critical exponent (f> is associated with the divergence of the specific heat at constant 
charge (Cq) at the critical point r + = r; L as, 



Cq ~ \Q - Qi 
Now from (32) and (36) we note that, 

1 



Cq 



which may be written as, 



Cq 



A 
1 



where we have used Equation (55). 
Comparison of (64) with (62) yields 

Critical exponent i/j: 



\Q-Q^ 2 
l 

2' 



(62) 

(63) 
(64) 



In order to calculate the critical exponent tp, which is related to the entropy of the third order 
LBI-AdS black hole, we Taylor expand the entropy (S(r+)) around the critical point r + = r, L . 
This gives, 

OS 



S(r+) = S(n) + 



(r+ - n) + 



(65) 



If we neglect all the higher order terms and use (19), (26) and (55), we can write (65) as, 



5-7T 3 

S(r + ) - S(r t ) = — (rl + 2ar* + a 2 ) 
Comparing (66) with the standard relation 



1/2 



\Q - Q, 



1 1/2 



5(r+) - S(n) ~ \Q - Q c 



(66) 



(67) 



we finally obtain ip = -. 

In the table below we write all the six critical exponents obtained from our analysis in a 
tabular form. For comparison we also give the critical exponents associated with some well 
known systems. 
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Critical exponents 


3rd order LBI-AdS black hole 


CrBr 3 * 


2D Ising model* 


van der Waals's system* 


a 


0.5 


0.05 










0.5 


0.368 


0.125 


0.5 


7 


0.5 


1.215 


1.7 


1.0 


5 


2.0 


4.28 


15 


3.0 


rP 


0.5 


0.60 






<t> 


0.5 


0.03 







(*: these are the non-mean field values.) 
(\: these values are taken from [36].) 



Thermodynamic scaling laws and static scaling hypothesis: 



The discussion on critical phenomena is far from complete unless we make a comment on the 
thermodynamic scaling laws. In standard thermodynamic systems the critical exponents are 
found to satisfy some relations among themselves. These relations are called thermodynamic 
scaling laws|261 127j. These scaling relations are given as, 



a + 2/3 + 7 
a + 13(5 + 1) 

<j> + 2V - \ 
13(5-1) 
(2-a)(5-l) 
l + (2-a)(5ip-l) 



2 
2 

1 

7 

7(1 + 5) 
(l-a)5 



(68) 



From the values of the critical exponents obtained in our analysis it is interesting to observe 
that all these scaling relations are indeed satisfied for the third order LBI-AdS black holes. 

We shall now explore the static scaling hypothesis [26j-[28j for the third order LBI-AdS black 
hole. Since we are working in the canonical framework, the thermodynamic potential of interest 
is the Helmhotz free energy, J-(T, Q) = M — TS, where the symbols have their usual meaning. 

Now the static scaling hypothesis states that, 

Close to the critical point the singular part of the Helmhotz free energy is a generalized 
homogeneous function of its variables. 

This asserts that there exist two parameters a e and ay such that 



•F(A ae e,A aT T) = AF(e,T) 



(69) 



for any arbitrary number A. 

In an attempt to find the values of the scaling parameters a e and ax, we shall now Taylor 
expand the Helmhotz free energy J-(T, Q) near the critical point r+ =r%. This may be written 
as, 



F(T,Q)=T(T,Q)\ 1 



+ 



df), 



\9Q 



+ 



(T - Ti) + 



(Q - Qi) + 

( d 2 F V 
\dTdQJ 



or 2 ) 
d 2 F 

dQ 2 



(T - Ti ] 



r+=ri 



(Q-Qi 



r+=ri 



(T-Ti)(Q-Qi) + 



(70) 



r + =ri 
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From (70) we can identify the second derivatives of J- as, 



and 



9t 2 ) q 



d 2 T 



-Cq 



(71) 



(72) 



dQ 2 ; , 

Note that, since both Cq and k~} diverge at the critical point, these derivatives can be 
justified as the singular parts of the free energy J 7 . 

Since in the theory of critical phenomena we are mainly interested in the singular part of the 
relevant thermodynamic quantities, we sort out the singular part of J-(T,Q) from (70), which 
may be written as, 



F, 



(T - Tif + 



1 



r+-- 



d 2 F 

dQ 2 



(Q - Qif 



-Cq 
2X; 



(T - Ti\ 



Qif 



(73) 



where the subscript V denotes the singular part of the free energy J 7 . 

Using (30), (36), (49) and (55) we may write the singular part of the Helmhotz free energy 
(J 7 ) as, 

F s = a^ 2 + r,T 3 / 2 (74) 



where 



and, 



n 



2 3/2 r 1 /2 



From (69) and (74) we observe that 



ax 



(75) 



(76) 



This is an interesting result in the sense that, in general, a t and ax are different for a 
generalized homogeneous function (GHF), but in this particular model of the black hole these two 
scaling parameters are indeed identical. With this result we can argue that, the Helmhotz free 
energy is an usual homogeneous function for the third order LBI-AdS black hole. Moreover, we 
can determine the critical exponents (a, f3, 7, 5, (ft, ift) once we calculate the scaling parameters. 
This is because these critical exponents are related to the scaling parameters as [2]3 |2"7] . 



a = 


2 

a e 


13 = 


1 — ax 

a e 


7 = 


2a T - 1 


a e 


5 = 


ax 


1 — ax 




eft = 


2a e - 1 




ax 


ift = 


1 — a e 


ax 



(77) 
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There are two other critical exponents associated with the behavior of the correlation function 
and correlation length of the system near the critical surface. We shall denote these two critical 
exponents as rj and v respectively. If G(f+) and £ are the correlation function and the correlation 
length respectively, we can relate rj and v with them as 

£~\T-Ti\~ v (78) 

and 

G{r + )~r 2 ~ n -\ (79) 
For the time being we shall assume that the two additional scaling relations [26] 

7 = v{2 - rj) and (2 - a) = vn (80) 

hold for the third order LBI-AdS black hole. Using these two relations ((80)) and the values of 
a and 7, the exponents v and n are found to be - and respectively. 

Although we have calculated n and v assuming the additional scaling relations to be valid, it 
is not proven yet that these scaling relations are indeed valid for the black holes. One may adapt 
different techniques to calculate 77 and i>, but till now no considerable amount of progress has 
been made in this direction. One may compute these two exponents directly from the correlation 
of scalar modes in the theory of gravitation [39J but the present theories of critical phenomena 
in black holes are far from complete. 



5 Conclusion 

In this paper we have analyzed the critical phenomena in higher curvature charged black holes 
in a canonical framework. For this purpose we have considered the third order Lovelock-Born- 
Infeld-AdS (LBI-AdS) black holes in a spherically symmetric space-time. We systematically 
derived the thermodynamic quantities for such black holes. We are able to show that some of 
the thermodynamic quantities (Cq, k^ 1 ) diverge at the critical points. From the nature of the 
plots we argue that there is a higher order phase transition in this black hole. Although the 
analytical estimation of the critical points is not possible due to the complexity of the relevant 
equations, we are able to determine the critical points numerically. However, all the critical 
exponents are calculated analytically near the critical points. Unlike the AdS black holes in the 
Einstein gravity, one interesting property of the higher curvature black holes is that the usual 
area law of entropy does not hold for these black holes. One might then expect that the critical 
exponents may differ form those for the AdS black holes in the Einstein gravity. But we find 
that all the critical exponents in the third order LBI-AdS black hole are indeed identical with 
those obtained in Einstein gravity [8] HI]. From this observation we may conclude that these 
black holes belong to the same universality class. Moreover, the critical exponents take the 
mean field values. It is to be noted that these black holes have distinct set of critical exponents 
which does not match with the critical exponents of any other known thermodynamic systems. 
Another point that must be stressed is that the static critical exponents are independent of the 
spatial dimensionality of the AdS space-time. This suggests the mean field behavior in black 
holes as thermodynamic systems and allows us to study the phase transition phenomena in the 
black holes. We have also discussed the static scaling laws and static scaling hypothesis. The 
static critical exponents are found to satisfy the static scaling laws near the critical points. We 
have checked the consistency of the static scaling hypothesis. Apart from this we note that the 
two scaling parameters have identical values. This allows us to conclude that the Helmhotz free 
energy is indeed a homogeneous function for this type of black hole. We have determined the 
two other critical exponents v and n associated with the correlation length (£) and correlation 
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function (G(r+)) near the critical surface assuming the validity of the additional scaling laws. 
The values of these two exponents are found to be - and respectively in the six spatial 
dimensions. Although the other six critical exponents are independent of the spatial dimension 
of the system, these two exponents are very much dimension dependent. 

In our analysis we have been able to resolve a number of vexing issues concerning the critical 
phenomena in third order LBI-AdS black holes. But there still remains some unsolved problems 
that encourage one to make further investigations into the system. First of all, we have made 
a qualitative argument about the nature of the phase transition in this black hole. One needs 
to go through detailed algebraic analysis in order to determine the true order of the phase 
transition[18|-[22j. Secondly, we have calculated the values of the exponents v and 77 assuming 
that the additional scaling relations hold for this black hole. But there is no evidence whether 
these two laws hold for the black hole[Hl These scaling relations may or may not hold for 
the black hole. The dimension dependence of these two exponents (y and rj) makes the issue 
highly nontrivial in higher dimensions. A further attempt to determine rj and v may be based on 
Ruppeiner's prescription [58J, where it is assumed that the absolute value of the thermodynamic 
scalar curvature (|7£|) is proportional to the correlation volume £ n : 

(81) 

where n is the spatial dimension of the black hole. Now if we can calculate 1Z using the standard 
method p2l 159 1 [60] . we can easily determine £ from (81). Evaluating £ around the critical point 
r + = T{ as before, we can determine v directly. It is then straight forward to calculate rj by using 
(80). This alternative approach, based on Ruppeiner's prescription, to determine v and rj needs 
high mathematical rigor and also the complexity in the determination of the scalar curvature 
(TV) in higher dimensions makes the issue even more challenging. 

Finally, it would be highly nontrivial if we aim to investigate the AdS / CFT duality as an 
alternative approach to make further insight into the theory of critical phenomena in these black 
holes. The renormalization group method may be another alternative way to describe the critical 
phenomena in these black holes. 
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